Dispersion relation of lipid membrane shape fluctuations by neutron spin-echo 

spectrometry 
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We have studied the mesoscopic shape fluctuations in aligned multilamellar stacks of DMPC 
bilayers using the neutron spin-echo technique. The corresponding in plane dispersion relation 
r~^(q||) at different temperatures in the gel (ripple, P/j') and the fluid (La) phase of this model 
system has been determined. Two relaxation processes, one at about 10ns and a second, slower 
process at about 100ns can be quantified. The dispersion relation in the fluid phase is fitted to 
a smectic hydrodynamic theory, with a correction for finite resolution. We extract values for, 
the bilayer bending rigidity k, the compressional modulus of the stacks B, and the effective sliding 
viscosity 773. The softening of a mode which can be associated with the formation of the ripple 
structure is observed close to the main phase transition. 

PACS numbers: 87.16.Dg, 83.85.Hf, 62.20.Dc 



Thermally excited shape fluctuations of membranes 
are readily observed in the light microscope, e.g., the 
flickering of lipid vesicles or biological membranes. On 
the nanometer scale, undulations of the lipid bilayer are 
experimentally much more difhcult to probe. Few spec- 
troscopic techniques offer the spatial and temporal res- 
olution needed for the validation of the theoretic work 
on collective bilayer motion, see e.g., Q, IMIMQilS' Up 
to now, thermal fluctuations of lipid membranes, in par- 
ticular in multilamellar stacks, have been mainly investi- 
gated by thermal diffuse scattering, i.e., x-ray lineshape 
analysis 0, Q- Elastic scattering has led to a detailed 
understanding of the static properties of thermal fluctu- 
ations in lipid membranes and the elasticity properties 
governing these fluctuations. According to linear smec- 
tic elasticity theory thermal fluctuations in the fluid 
phase of the membrane are governed by the free energy 
functional (Hamiltonian) 
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where k denotes the bilayer bending rigidity, A the area 
in the a;y-plane, N the number of bilayers, and Un the de- 
viation from the average position n d of the n-th bilayer, 
d is the lamellar spacing. B and K = n/d are elastic coef- 
ficients, governing the compressional and bending modes 
of the smectic phase, respectively. A fundamental length 
scale in these systems is given by the smectic penetration 
length A = y^K/B. Aligned lipid bilayers allow a sepa- 
rate determination of both parameters K and B yjJ, [llj . 
However, a full understanding of the collective dynamics 
should include experimental determination of relaxation 
rates, the relevant transport coefficients (viscosities), and 
in particular the characteristic dispersion relation of the 
relevant modes. 



In this letter we present an experimental dispersion re- 
lation for aligned multilamellar lipid membranes T~^{q\^) 
as a function of lateral momentum transfer (in the 
plane of the bilayers), measured by neutron spin-echo 
spectrometry (NSE). From analysis of the dispersion rela- 
tion, the effective sliding viscosity, 773 , as well as the static 
properties k and A are obtained. A softening of a mode 
which can be associated with the formation of the ripple 
structure is observed close to the main phase transition. 
In the (7|| range probed, the bilayer displacement u„ can 
be assumed to be a continuous variable. It is therefore 
not sensitive to the discrete molecular structure, in con- 
trast to the dynamics on the molecular length scale (acyl 
chain distance) which can be probed by inelastic neutron 
scattering using three-axes spectrometry We have 

selected NSE for this study since the undulation modes 
at high g|| are too fast to be accessed by x-ray photon 
correlation spectroscopy (XPCS) and the lateral length 
scales are to small to be resolved by dynamic light scat- 
tering (DLS) [3. XPCS has successfully been used for 
investigations of free-standing thermotropic liquid crys- 
tals, which exhibit much slower acoustic modes, corre- 
sponding to center of mass movement of the entire film 
|l4j . NSE has been used before to study phospholipid 
membranes, |l5lll6l |. However, a dispersion relation over 
a wide range in g|| in the fluid multilamellar state has 
not been achieved so far, in part due to flux and techni- 
cal limitations in the past. 

NSE directly measures the intermediate scattering 
function S(q;t), which is related to the density distri- 
bution p(r, t) in the sample by 

Siq, t)= [ d^Re-^qR [ d'r {p{r, 0)p(r + R, t)) . (2) 



The spin-echo experiments were carried out at the INll 
and IN15 spectrometers, situated at the cold source of 
the high flux reactor of the Institut Laue-Langevin (ILL) 
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FIG. 1: (Color online), (a) Sketch of 
the scattering geometry, (b) Rock- 
ing scan at T=30°C through the 
first Bragg peaks and along the dif- 
fuse Bragg sheet in log-log scaling. 
The arrows exemplary mark posi- 
tions of inelastic scans in the fluid 
phase. (c) Inset shows a reflec- 
tivity curve as measured on INll 
with relaxed momentum resolution 
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FIG. 2: (Color online). Intermediate scattering function S(q,t)/S(q,0) for selected qj] values in the interval 
0.002A~^ <q|l <0.8A~^ measured on (a) INll at T=22°C, just above the temperature of the main phase transition in 
deuterated DMPC -d54, (b) on IN15 at 30°C, in the fluid phase of the phospholipid bilayers, and (c) on IN15, at 19°C in the 
gel phase. Two relaxation processes, one at about 10ns and a second slower process at about 100ns are clearly distinguished. 
Solid lines in the Figures represent least-square fits to two exponential decays. 



in Grenoble, France. Wavelength bands centered at 
A=7.4A and A=14A with AA/A ~ 0.15 (FWHM), re- 
spectively, have been set by a velocity selector. Partially 
(acyl chain) deuterated DMPC-d54 (1,2-dimyristoyl-sn- 
glycero-3-phoshatidylcholine) was obtained from Avanti 
Polar Lipids. Highly oriented multi lamellar membrane 
stacks of several thousands of lipid bilayers were pre- 
pared by spreading lipid solution of typically 25mg/ml 
lipid in trifluoroethylene/chloroform (1:1) on 2" silicon 
wafers, followed by subsequent drying in vacuum and hy- 
dration from D2O vapor |l7| . resulting in a structure of 
smectic A symmetry. Twenty such wafers separated by 
small air gaps were combined and aligned with respect 
to each other to create a "sandwich sample" consisting 
of several thousands of highly oriented lipid bilayers (to- 
tal mosaicity about 0.5°), with a total mass of about 
400mg of deuterated DMPC. During the experiment, the 
samples were kept in a closed temperature and humidity 
controlled Aluminum chamber, and were hydrated from 



the vapor phase. 

The layer fluctuations of lowest energy arc the un- 
dulation modes, i.e., highly correlated (conformal) layer 
displacements which keep the inter-layer distances ap- 
proximately constant. This conformality explains the 
strong modulation of the diffuse scattering in q^, i.e., 
the Bragg sheet structure 0, corresponding to the 
diffuse tails (at > 0) of the specular Bragg peaks 
(qz = n 2-K/d,q\\ = 0). We have now measured spin- 
echo curves on the first Bragg sheet as a function oi q\\, 
i.e.; at constant q^ = 27r/(i. A sketch of the scattering ge- 
ometry is shown in Fig.n(a). The in plane component of 
the scattering vector is calculated to q||=q2 tan(tt;), with 
q^ w O.llA^^ for the first reflectivity Bragg peak, lo is 
the rocking angle, i.e., the sample rotation with respect 
to specular Bragg angle. Figure n](b) depicts a rocking 
curve in the fluid phase at T=30°C on the first reflectiv- 
ity Bragg peak where the sharp Bragg component and 
the diffuse Bragg sheet can be well distinguished. In the 
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inset, (c), a reflectivity curve is plotted. Peaks are dis- 
tinctly broadened by the wavelength distribution. 

Figure El shows the intermediate scattering function 
S(q||,t) for selected q|| values for spin-echo times of 
0.001ns<t<20ns for INll and 0.01ns<t<200ns for IN15. 
Data have been taken at three different temperatures, at 
19°C, in the gel (ripple, Pp') phase of the phospholipid 
bilayers, at 22°C, just above the temperature of the main 
transition in deuterated DMPC-d54 (at « 21.5°C), 
and at 30° C, far in the fluid Lq phase of the mem- 
branes and above the regime of so-called anomalous 
swelling. The corresponding lamellar d spacings were 
d=56k, 60 A and 54A (gel, 22°C and fluid), respectively. 
Two relaxation processes, one at about 10ns and a sec- 
ond, slower process at about 100ns are clearly distin- 
guished. Solid lines in the Figure represent least-square 
fits of the data to two stretched exponential decays, 
5(g||,t)/5((7||,0) = {A, - A2)cxph {t/n{qii)f']+y2 + 

(^2 -2/2)exp[- (i/T2(i?||)) ]. The fitting results for /3i^2 
were in the range of 0.96 < /3i_2 < 1 for the temperatures 
T=22°C and 30°C, corresponding essentially to single 
exponential relaxations. Contrarily, the gel phase results 
(19°C) give values of /3i^2 ~ 1.75, i.e., compressed expo- 
nentials. Structural inhomogeneities and heterogeneous 
interactions would lead to a local relaxation dynamics 
and to stretched (/3 < 1) exponentials. A compressed ex- 
ponential decay is incompatible with a diffusive, fluid like 
motion of the particles and might therefore be an intrinsic 
property of the gel state. The relaxation rates t^^ and 
T2~^ in the gel and the fluid phase are depicted in Fig.|2| 
(a) and (b), after compilation of all measured q|| values 
on the two spectrometers. Also shown is the fast relax- 
ation process for T=22°C. Both relaxation branches are 
clearly dispersive. The fast process shows a q^| increase 
at small q|| values and a bend at about q|| «0.015A^^. 
The dispersion in the gel phase and close to the phase 
transition in Fig. 13 (b) appear to be more pronounced as 
compared to 30°C dispersion. As a remarkable feature, 
a soft mode appears in the T=22°C dispersion, indicat- 
ing a significant softening of the bilayer at a well defined 
wave number. The slow branches at T=19°C and 30°C 
also show increasing relaxation rates with increasing q|| 
values, but with a distinct non-polynomial behavior. 

The dispersion relation of the fluid Lq phase lends itself 
to a quantitative comparison with theory. As an effect 
of the stacking of the bilayers and the substrate, undula- 
tion modes in lamellar phases decay with relaxation rates 
T^^((j'll) = K/{ri^d) g^, while fluctuations of a free, non 
supported, bilayer would decay with °^ ^j] S 
Generally, undulation modes are probed at q^ values of 
qz = 27: /d. If the scattering is probed at flnite com- 
ponents 5qz '■= {qz — 27r/(i) or measured with a relaxed 
qz resolution, baroclinic modes rather than pure undula- 
tions are probed. Dynamic light scattering experiments 
can be described according to bulk smectic elasticity the- 
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FIG. 3: (Color online), (a) Dispersion relations at T=30°C. 
The solid line is a fit to Eq. (b) Dispersion relations in the 
gel (19°C) and in the fluid phase (22°C). A pronounced soft 
mode is observed at qo «0.015A~^ at 22°C (dotted vertical 
line). Solid lines in (b) are guides to the eye. 



ory and multi lamellar fluctuations as a function of both 
symmetry axis {q\\,qz) by the dispersion relation |l8l |: 

, Bql + (K/d)q 



"(911 



The transport coefficient ji 
12770) with ?7o the viscosity of 



^11 1)39' 

was estimated to fj, = d^ 

the the solvent (water) [ijj . Note that higher order elas- 
tic terms have been neglected. The neutron scattering 
experiments are carried out in the first Brillouin zone so 
that we have to replace q^ by 6qz , the deviation from the 
first zone center. However, finite size or resolution ef- 
fects lead to a broadening of the diffuse Bragg sheet and 
therefore set a lower cut-off in the reachable Sqz ■ Follow- 
ing the idea of Ribotta the above equation is then 
replaced by 



_i _ n/d 2<ll\ 



{7r/{AD)Y 



(3) 



In our case, the qz width of the diffuse scattering is not 
defined the finite size D of the lipid film, but comes from 
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the rather broad instrumental resolution, which is given 
by the wavelength spread. From the measured width 
of the first reflectivity Bragg peak in Fig. ^ (c) of 
Aq^ ~ 0.015A-1 (HWHM, Lorentzian fit), we obtain an 
effective finite-size cutoff-length D = tt/ Aqz ~ 212A that 
we used to fit the measured fluid dispersion relation to 
Eq. see solid line in Fig. 13 (a). The component of 

the scattering vector was neglected, q = + gj^ ~ 9| | , 

because it may average out for symmetric broadening. 
Alternatively, the results can be compared to the theory 
given in |2| for thin film samples, including the bound- 
ary effects and surface tension. However, this theory 
gives dispersion relations for each of the {N — 1) eigen- 
modcs, which makes a comparison difficult if not impos- 
sible without proper weighting of the modes. 

Using Eq. JSJ, the following results are obtained for 
the three free parameters: k = 14.8 ± Sk^T, A ~ 
10.3 ± 2.3A, rys = 0.016 ± O.OOOGPa s. B is calculated 
to B=1.08 lO^J/m^ (d=54A). Note that the value for A 
compares quite well with the value obtained by a very 
different approach but for a similar swelling state [20|. 
Experiments at different osmotic pressures point to a 
distinct effect of the swelling state to the compressional 
modulus B and A, consequently 0. The bending modu- 
lus compares quite well to results from Molecular Dynam- 
ics although the simulations do not reach the small q| | 
values that we probe experimentally, k has large errors, 
but lies in the middle of the even larger range of liter- 
ature values 113. q|| values of q|| >0.05A~^ in Fig. |31 
(a) and (b) show declining relaxation rates which devi- 
ate from the theoretic curve in Fig. |21 (a). The diffuse 
Bragg sheet was slightly bend for the high q| | values and 
the corresponding diffuse scattering was likely dominated 
by defect and no longer purely thermal diffuse scattering 
with consequences for the precision in the determination 
of K and 773. The dynamics at low q\\ values is governed by 
the interplay of viscosity and inertia, i.e., r^^ cx Tj^/pq^^, 
while at higher q|| values, the undulation dynamics be- 
comes predominant. The difficulty in quantifying relax- 
ation times at higher qj[ values thus limits the precision 
in the determination of k distinctly. This constraint can 
be overcome in future studies with optimized set-ups and 
sample preparation. Note that conceptually, fluctuations 
in the gel phase, should not be described by the disper- 
sion relations of fluid smectic A phases. The solid lines 
in Fig. 13 (b) are therefore guides to the eye. 

In summary, NSE provides unique access to collective 
dynamics of shape fluctuations in solid supported mul- 
tilamellar lipid membranes. The dispersion relation of 
the fast branch with relaxation rates between 1 and lOns 
can be attributed mainly to undulation dynamics with 
the expected mixing of baroclinic modes resulting from 
a correction for finite resolution. The resulting effective 
sliding viscosity of the membrane system 773 was found to 
be 16 times higher than that of water. The slow disper- 



sion branch with relaxation rates of about 100ns may be 
attributed to a surface relaxation mode and will be 
addressed in a forthcoming publication. Furthermore a 
quite localized deviation from the undulation branch was 
observed at T=22°C just above the main phase transi- 
tion, indicating the softening of a well defined wave num- 
ber of qo ~ O.OISA"^. From this result, we speculate that 
the well known softening of phospholipid membrane upon 
approaching the main phase transition temperature from 
the fluid phase, i.e., the regime of "critical swelling" or 
" anomalous swelling" , occurs on a well defined length 
scale, here at 27: /q^ « 420 A. We have additionally mea- 
sured the elastic scattering in the ripple phase and from 
the satellite peaks of the diffuse Bragg sheet we deter- 
mine a ripple periodicity of dj- w 130A, distinctly smaller 
than the length scale of the soft mode observed here. 
However, dr grows significantly at the transition, as we 
have measured by atomic force microscopy. Therefore, 
we speculate that the soft mode in the fluid phase is 
linked to the formation of the ripple structure in the P^i 
phase, but this remains to be investigated in more detail. 
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